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Abstract 



A theory based on the thermodynamic Gibbs-Thomson relation is pre- 
sented which provides the framework for understanding the time evolution of 
isolated nanoscale features (i.e., islands and pits) on surfaces. Two limiting 
cases are predicted, in which either diffusion or interface transfer is the limit- 
ing process. These cases correspond to similar regimes considered in previous 
works addressing the Ostwald ripening of ensembles of features. A third pos- 
sible limiting case is noted for the special geometry of "stacked" islands. In 
these limiting cases, isolated features are predicted to decay in size with a 
power law scaling in time: ^ oc (to — t)"', where A is the area of the feature, 
to is the time at which the feature disappears, and n = 2/3 or 1. The con- 
stant of proportionality is related to parameters describing both the kinetic 
and equilibrium properties of the surface. A continuous time Monte Carlo 
simulation is used to test the application of this theory to generic surfaces 
with atomic scale features. A new method is described to obtain macroscopic 
kinetic parameters describing interfaces in such simulations. Simulation and 
analytic theory are compared directly, using measurements of the simulation 
to determine the constants of the analytic theory. Agreement between the 
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two is very good over a range of surface parameters, suggesting that the an- 
alytic theory properly captures the necessary physics. It is anticipated that 
the simulation will be useful in modeling complex surface geometries often 
seen in experiments on physical surfaces, for which application of the analytic 
model is not straightforward. 
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I. INTRODUCTION 



Real surfaces of crystalline materials below their roughening temperature, even of crystals 
with essentially perfect ordering, are generally not in equilibrium. More typically they have 
features such as irregular step edges, excess individual adatoms or vacancies, or clusters 
of adatoms or vacancies which we call islands and pits. Since all these features have an 
associated free energy cost but are not required by the macroscopic orientation of the surface, 
their existence indicates that the surface is not at a free energy minimum. Such surfaces 
will therefore tend to relax towards equilibrium. 

The way in which this relaxation takes place has been a subject of interest for many years. 
In processes which tend to produce nonequilibrium features, such as epitaxial growth or 
surface sputtering, the relaxation may be important in determining the time evolution of the 
surface both during and subsequent to deposition or irradiation. Surface relaxation will also 
be of particular interest in the future for the field of nanofabrication. As artificially created 
structures approach the nanometer scale, it becomes important to understand whether these 
structures, once made, will be relatively stable. Stability of such features is crucial to success 
for electronic interconnects or information recording applications. 

One of the most studied situations is a surface populated by an ensemble of islands of 
various sizes, as might be found, for instance, in a growth experiment after the initial stage 
of island nucleation. This is a special case of a more general situation in which small clusters 
of a condensed phase of material, which are of dimensionality D^., exist in an environment of 
a vapor phase, which is of dimensionality D^. Relaxation towards equilibrium in many such 
systems occurs through the diffusion of atoms between the clusters, in a behavior known 
as Ostwald ripening (also called coarsening). Greenwoodi and Lifshitz and SlyozovHI first 
considered the case Dc = D,„ = 3 under the simplifying assumptions that the inter-cluster 
mass flow is limited by the diffusion process, and that the total amount of the condensed 
phase is small. WagnerS extended the theory to include the possibility that the transfer of 
adatoms between the phases could be a limiting process. Further extensions by ChakravertyB 
and others!^ applied the theory to the case of clusters on surfaces (i.e., Dy = 2, Dc = 2 or 
3). References and |TT] include reviews of this material. 

All of these treatments and many experimental investigationsllUHl have concentrated on 
quantities, such as average cluster radius and cluster density, which are ensemble properties. 
This is because the clusters involved are usually very small, so that their properties are most 
readily studied by techniques, such as diffraction, which sample large numbers of clusters. 
With the advent of atomic scale microscopies such as scanning tunneling microscopy (STM), 
it has become practical to study the behavior of individual surface features, either within an 
ensemble or in isolation. By observing the behavior of individual features, the fundamental 
mechanisms responsible for stability and decay can be studied without the complications 
which result from the interactions between the clusters of an ensemble system. In addition, 
the behavior of an individual feature may be the quantity of interest in certain applications, 
as in nanofabrication. Some studies have focused on the evolution of individual features, 
either experimentallyil^0 or theoretically.^^ For instance, we have found that, on the 
Au(lll) surface in air, isolated islands decay in size such that the island area depends 
roughly linearly on time.H A similar study of islands on the Ag(lll) surface in ultra high 
vacuumil observed island areas to be proportional to time raised to an average exponent of 
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0.54. 

In this work we present an extension of the ensemble theories to consider the decay 
of isolated individual two dimensional islands on surfaces {Dy = = 2). In Section |I| 
an analytic continuum model is presented which accounts for three possible rate limiting 
steps in island decay. This model provides a connection between the observed behavior of 
such features and the parameters describing atomic processes on the surface. Section |T| 
discusses limiting cases in this model which are relevant to physical systems. Specifically, 
we find "interface-transfer limited" and "diffusion limited" cases which are direct corollaries 
of similar limits of Ostwald ripening.i0 Specific surface geometries lead to cases which are 
unique to individual islands. Islands completely surrounded by a step down or step up edge 
can lead to the "outer-boundary limited" and "island in a pit" cases. Villain has considered 
the first of these geometries in Ref. |2^, using assumptions which lead to the diffusion limited 
case. 

Because the evolution of a single surface feature involves a relatively small area of the 
surface, it is feasible to perform atomic scale simulation studies to test the analytic model. 
Section presents such a computer simulation study of island decay for a generic surface 
with one particular geometry. In Section we discuss how to obtain the parameters for the 
continuum model which are appropriate to describing the atomistic simulation. Using these 
parameters, we find in Section |V^ that comparison between the simulations and the analytic 
model is very good, showing that the analytic model captures the key physics for describing 
feature evolution on a generic surface. We anticipate that this simulation technique will 
provide an additional tool for analyzing the observed behavior of physical surfaces, because 
it can be used to model complex geometries which are often observed, e.g., in STM studies, 
but which are not easily treated with the analytic model. 



II. THEORY FOR DECAY OF ISOLATED SURFACE FEATURES 
A. The driving force for surface mass flow 

At the heart of each of the treatments cited above is the Gibbs-Thomson (GT) relation. 
The pressure of a vapor which is in equilibrium with its condensed phase depends on the 
curvature of the interface between the phases. For a convex condensed phase the pressure 
is enhanced relative to the pressure in equilibrium with a flat boundary. The GT relation 
expresses this dependence. The increased vapor pressure associated with a curved boundary 
is the driving force behind the mass transfer between different clusters or between a cluster 
and a nearly straight step edge. 

For a two dimensional, incompressible condensed phase, the boundary of which has a 
radius of curvature r, in contact with a two dimensional, noninteracting vapor, the GT 
relation isS 

where p'^'^ is the pressure at which the vapor will be in equilibrium with the condensed 
phase, is the equilibrium vapor pressure for a straight boundary (r = oo), 7 is the edge 
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free energy (which is assumed to be isotropic and independent of curvature), and Q is the 
area occupied by one atom in the condensed phase. When working with a lattice, it is more 
natural to consider the vapor concentration p, measured in units of where a is the lattice 
constant. Because Eq. (|lD assumes that the vapor is noninteracting, the pressure will be 
proportional to the vapor concentration. Thus we can substitute p for p in Eq. (|1]). 

Note that the same relation continues to hold for a condensed phase with a concave 
boundary. In that case, the boundary is described as having a negative radius of curvature. 
For example, the adatom density inside a monolayer deep pit will be reduced relative to the 
pressure in equilibrium with a fiat boundary. 

If the vapor is interacting, as is the case for adatoms on a surface, we have reported 
elsewhere^ that at vapor concentrations as low as p^"^ ~ 10"'^ there can be significant 
corrections to this equation. However, by including a correction factor in the exponent 
an equation of the same functional form as Eq. (0) continues to satisfactorily describe the 
Gibbs-Thomson vapor enhancement in our system, as will be discussed further in Sec. |IV|. 

The Gibbs-Thomson relation is sufficient to anticipate the general behavior of a system. 
Features with the highest curvature will have the largest equilibrium adatom vapor concen- 
trations. On a surface with islands of different sizes as well as step edges, a concentration 
gradient will be established whereby adatoms diffuse away from high curvature features and 
toward low curvature features. For instance, large islands will grow as small islands shrink 
and disappear. 

We now consider in more detail the specific case of an isolated circular island on a 
surface with an adatom concentration below the island's equilibrium value. We shall use r 
to denote the radius of the island, and r to denote the radial distance from the center of 
the island. The GT equation describes the equilibrium properties, but in order to calculate 
the decay rate of an isolated island we must consider the kinetics of the surface. In the 
following sections, we will describe the three processes involved in the flow of atoms from 
the island: interface transfer ,&i diffusion,!! and incorporation at the outer boundary.! These 
three processes are illustrated in Fig. |l|. We make the assumption throughout this treatment 
that these processes are slow on the time scale of local equilibration in the vapor, so that 
the system may be described by a steady state solution at any point in time. 



As the island decays, there must be a net detachment of atoms from the island. If the 
deviation of the density of atoms just outside the island p{r = r) from the equilibrium density 
given by the GT relation is sufficiently small, then the current of atoms leaving the island 
Jj will depend linearly on that deviation. Further, if the curvature of the interface does not 
significantly influence the manner in which atoms leave the island then the current will be 
proportional to the interface length, so that 



where we introduce the coefficient Ki which is assumed to be independent of r. The subscript 
'i' is chosen to identify the net detachment current as the interface transfer current of atoms 
out of the island. 



B. The interface transfer process 




(2) 



5 



Further insight into the coefficient Ki can be obtained from a detailed balance argument. 
In the remainder of this subsection we explore how this coefficient arises from the microscopic 
dynamics, thus clarifying the assumptions made above and forming a basis for predicting 



the behavior of physical systems in Sec. |III B| . In order for an atom to attach to or detach 



from the island, it must overcome some energy barrier at the island-vapor interface (which 
depends on the particular microscopic atomic move). Let Jd be the detachment current 
density at the interface, i.e., the rate per unit length of perimeter at which atoms detach 
from the island edge. Similarly, let Ja denote the interface attachment current density. Both 
of these will, in general, depend on the interface curvature and the vapor concentration. If 
we assume that both p(r) and [p^^ — p{r)] are ^ 1 [an assumption more stringent than 
required for Eq. (||), but frequently satisfied] then we may proceed by Taylor expanding the 
current densities in p{r), obtaining 

Ur,p{r)) = J'Ar) + ^^(r)p(r) + \j'l{r)[p{r)f + • • • , (3a) 

Ja{r, P{r)) = mp{r) + \j:{r) [p{r)f + ■ ■ ■ , (3b) 

where J^if) is the detachment current density with no vapor, and the prime indicates 
partial derivative with respect to p evaluated at p = 0. Note that there is no constant 
term in Eq. (^) because there can be no attachment unless an adatom vapor is present. 
If the coefficients are well behaved we can drop terms of higher than linear order. Under 
equilibrium conditions, when p(r) = p^^, the rate of detachment and attachment are equal, 

j'Ar) + ^^(r)p^.^ = J:(r)p^,^ . (4) 
The total current /j is determined by the current densities, which simplifies through Eq. (^, 

li = 27rr {Jd{r,p) - Ja{r, p)) 

= 2nr Mr) - J^(r)) (p^,^ - p(r)) . (5) 




R 



FIG. 1. A schematic representation of the adatom vapor density near a small island of radius r 
on a surface. The equilibrium vapor densities for the island (pf^) and for the mean field environment 
at radius R (pr) are indicated by dashed lines. The actual density profile p{r) lies between these 
limits, determined by the relative rates of the three processes of decay indicated by arrows. 
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Comparison with Eq. (Q) yields an expression for Ki, 



K^ = fair) - J'Ar) ^ J'air) , (6) 

where we have used the fact that in general we expect the detachment rate Jd to have 
little dependence on the vapor concentration, i.e., | -C \ Ja\- Equation gives a clearer 
meaning to the assumption that Ki is independent of curvature; it is equivalent to the 
assumption that J'a{r) (qualitatively the rate at which an atom will attach to the island 
if it is already one step away, i.e., a rate per adatom density) is independent of r. This 
assumption is plausible; for instance, if we have an ideal gas in an enclosure, the number of 
gas atoms hitting a given area of the enclosure wall is independent of its curvature. However, 
this assumption should be justified for any particular system under study. 



C. The diffusion process 



Once the atoms have detached from a decaying island, they must diffuse away across the 
terrace. Under our assumption of quasi-steady state, the rate at which they diffuse away 
is governed by the two-dimensional (2D) time-independent radial diffusion equation. The 
general solution of this equation diverges logarithmically at large distances. It is therefore 
necessary to specify a boundary condition far from the island. It is assumed that the island is 
sufficiently isolated from other surface features so that a mean field approach is valid. Rather 
than specify the island's environment in detail, we apply a generic boundary condition that 
the vapor concentration is some p{r = R) at a radius R> r, where R may vary in time. We 
will consider some particular examples below. The current of adatoms is then given by the 
usual solution to the steady state 2D radial diffusion equation: 

= 2nKD(fi{r) - p(R)) (7a) 



where we introduce the coefficient Kd- D is the collective (or chemical) diffusion coefficient 
of the vapor on a step-free terrace. The subscript 'D' is chosen to denote diffusion. 



D. The process of incorporation at the outer boundary 



Finally, the adatoms must leave the vapor by crossing the outer boundary, thus becoming 
incorporated into the mean field. Since we have left unspecified the nature of the outer 
boundary, this process may limit the rate at which atoms can leave. We may reason in a 
manner similar to the case of interface transfer, as discussed above. There will be some 
concentration such that if p{R) = pr then the net current across the outer boundary 
vanishes. If p{R) does not vary much from pR, then the current across the outer boundary 
will be proportional to their difference. 
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h = 27iRKb{piR) - pr) 



(8) 



This equation introduces the rate coefficient K^, where the subscript '6' is chosen to denote 
the outer boundary. Note the similarity between Eqs. (|) and (0). 



E. The total mass flow and island decay rate 



As commented at the end of Sec. [11 A] , we assume that the system is in a steady state. 
Therefore the three atom currents, which are in series, must be equaL Equating Eqs. (^, 
([7a|), and (||) gives the following expressions for the concentrations of adatoms at the bound- 
aries and for the net detachment current in terms of equilibrium concentration values: 



2iC(r) - p^) 



where 



C{r) 



exp 



7^ 
rkT 



PR_ 

eq 
Po2 



Kir 



1 



(10a) 
(10b) 

(11) 



In the Eq. (|10b|) , the GT relation has been substituted for p'^^. In Eq. (|TT]) the coefficients 
of the individual current equations add together as conductances. Thus C(r) is the overall 
conductance relating the adatom current to the driving force provided by an imbalance in 
equilibrium densities. Note that p{r) and p{R) are intermediate between the equilibrium 
concentrations pl^ and pr (see Fig. |I]). In this way a driving force is present for all three 
processes in the mass flow away from the island. 



III. LIMITING CASES OF THE ANALYTIC MODEL 

A. Simplifying mathematical limits 

The form of Eqs. ([I0|) and ([TT|) is quite complex, especially in their dependence on 
r. In order to better understand them it is useful to consider simplifying mathematical 
limits. Using physically reasonable approximations, we will find limits in which the decaying 
island size obeys a simple power law. That is, in these limits the area of the island obeys 
A oc (to — t)" where A is the area of the island and to is the time at which the island 
disappears. It will be seen that the exponent n may take the values of 2/3 or 1. 
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Two approximations must be satisfied in order to achieve either of the power law hmits. 
The first approximation is that (i) the equihbrium adatom density at the outer boundary 
is close to the equilibrium density for a straight step edge, pn ~ p^. We expect this to be 
valid if the island edge has a much higher curvature than other nearby surface features, so 
that Ip^'^ — pJJ^I ^ |p/j — p^|. We will thus label this as the "i? ^ r" approximation. 

The second approximation required is that (ii) the argument of the exponential in 
Eq. ( |10bD must be much less than one, so that the exponential may be approximated by the 



constant and linear terms from its Taylor expansion. This approximation we will label as the 
"GT expansion". It is also frequently physically reasonable. For instance, for a surface with 
7 = 0.117eV/a (see Sec. fV[) at room temperature, approximation (ii) is good to within 10% 
for islands of size r ^ 23 a. This approximation is not always valid, however; for instance, 
it seems to fail for island decay observed by Morgenstern et a/.0 
Using these first two approximations, Eq. (|10b|) may be rewritten 

/^27rC(r)p«^^. (12) 

In order to get a power law decay of the island size, a third required approximation is that 
(iii) C(r) must be constant or proportional to a power of r. This can be realized in two 
ways. (Physical realizations of the following approximations will be discussed in Sec. [ill lj| .) 



In one limit (iii-a) C(r) may be approximated to be a constant Cq. This can occur, for 
instance, if the third term of Eq. ([TT|) is much larger than the first two, i.e., K}jR <^ min(i^'jr, 
Kj:,), and if R is constant in time. This we call the "constant C" approximation. The adatom 
current gives the time rate of change for the island area, so that in this case we have 

This may be directly integrated to obtain 

r = {^Co^^y\to - ty/\ (13b) 

A oc (to - t)" , n = 2/3 . (13c) 

This will be referred to as the "ra = 2/3 mathematical limit". 

In the other limit (iii-b) of Eq. |ll] which can lead to power law behavior, C{r) is ap- 
proximately proportional to r. This we call the "linear C" approximation. This can occur, 
for instance, if the first term of Eq. { ^A\j is much larger than the last two. Expressing the r 
dependence of C(r) explicitly with C(r) = Cir, analogously to Eqs. (|13D we have 

2vrr| = ^ = -m = -2vrf^Cip-^ , (14a) 
r=(2C,P^^Y\to-tY/\ (14b) 
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A oc (to - t)" , n = l 



(14c) 



In this "n = 1 mathematical hmit" , the island area scales linearly with time. 

Note that in both of these limits, the velocity of the interface dr/dt diverges as the island 
approaches the vanishing point at time to- If £1 high enough velocity is reached, then the 
assumption that the system is in a quasi-steady state will fail. In practice, the velocity is 
limited by the fact that the island cannot be smaller than one atom, so that the breakdown 
of quasi-steady state might not occur. 

The formalism presented above for describing the decay of islands via adatom diffusion 
may also be applied to the decay of pits in the surface (vacancy clusters) via vacancy diffu- 
sion. These pits may be the source of vacancies which diffuse in the surface and recombine 
at step edges. In principle, both adatom and vacancy motion can contribute to mass flow 
on the surface. However, the rate constants may be very different for adatom and vacancy 
motion, so that one may dominate the surface evolution. 



B. Particular examples for physical systems 

In this section we will discuss limiting behaviors associated with each of the three pro- 
cesses involved in the mass flow from a decaying island, i.e. interface transfer, diffusion, and 
outer boundary incorporation (see Sees. [11 B| , [11 C| , and [11 D[ respectively). The special case 
of an island located within a pit will also be discussed. 

The treatment of Sec. ^ bears similarities to treatments concerned with Ostwald ripening 
in ensembles of clusters on surfaces.li0 In those treatments, a central quantity of interest is 
the critical radius Vc, defined such that islands of size r < Tc shrink while larger islands grow. 
Under approximations analogous to those in Sec. [Ill A[, along with certain assumptions about 



the ensemble size distribution, the critical radius is found to increase in time according to 
a power law.0 Indeed, it can be shown that the resulting exponents are identical to those 
found above for an individual shrinking island.ii Note, however, that the critical radius of 
an ensemble grows, while the physical radius of an isolated island shrinks. 

The physical characteristics which lead to power law behavior in Ostwald ripeningilisl will 
also lead to corresponding behavior for isolated islands. Consider a configuration in which 
an island is not surrounded by any clearly defined outer boundary. For instance, it might 
be on a terrace bounded on two sides by relatively straight steps. We assume that the outer 
boundary incorporation process described in Sec. [II D does not significantly affect the island 



decay, so that the third term in Eq. (^) is negligible. If the curvature of the island edge is 
higher than that of other surface features, i.e., if the island is small, then the approximation 
(i)(i? ^ r) remains valid. The island must be larger than than the length scale set by the 
step free energy, so that the approximation (ii)(GT expansion) is also valid. The choice 
of a mean field outer boundary R is problematic (it still enters the problem through Kd). 
If it is very difficult for atoms to cross from the condensed phase to the vapor phase and 
vice versa at the island edge, then the first term of Eq. ([TI[) dominates the others. In this 
way C{r) K^r, so that the outer boundary is totally removed from the problem. The 
approximation (iii-b) (linear C) and Eqs. (|l^) apply. This is called the "interface-transfer 
limited" case, and falls into the n = 1 mathematical limit. 
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If, on the other hand, the isolated island is in a system where diffusion is very slow, 
then we have the "diffusion limited" case, for which the second term of Eq. (pUD dominates. 
Since C(r) ^ depends on the outer boundary through Eq. ([75D, further assumptions 
must be made. The classic treatment of this regime! assumes that the weak logarithmic 
dependence of on r and R may be neglected. This is equivalent to assuming that there 
is a vapor phase screening length ^ = {X — l)r. This screening length leads to an outer 
boundary that goes as R = Ar, so that C{r) ^ Kd = D/ln{X) is constant. Thus the 
approximation (iii-a) (constant C) and Eqs. ([T3|) apply, and the diffusion limited case falls 
into the n = 2/3 mathematical limit. It should be noted, however, that the screening length 
hypothesis is not physically self-consistent, so that this approach should not be expected to 
yield quantitatively correct results. In particular, a screening length can only arise through 
the interaction of many features (islands, pits, and step edges), in which case it turns out 
that the screening length is too long to properly treat the island as isolated.i'ii 

A cross over between the interface-transfer limited and diffusion limited cases can be 
defined by the condition that the first two terms of Eq. ([IT]) are equal. At sufficiently large 
radii KiV will be larger than Kd, and the mass flow must be in the diffusion limited regime. 
Only when the island reaches a sufficiently small size can the interface-transfer limited regime 
be observed. Physically, a shorter island perimeter provides fewer detachment sites, making 
interface transfer more difficult. We can define a transition radius at which this cross over 
will occur by equating KiV and K^, which gives [see Eqs. (^ and ([TbD] 



Since ln(A) is a number of order one, we see that in order for to be significant, must be 
significantly larger than D. That is, there must be some impediment to an atom attaching 
to an island relative to diffusion. On a clean metal surface this is unlikely, as attachment 
energy barriers are generally found to be similar to or lower than diffusion barriers. @@ For 
adsorbate covered surfaces, however, the adsorbate will often tend to bind to the step edge 
and may therefore present a steric impediment to attachment. 

For individual islands, other cases arise which have no analog in the ensemble treat- 
ments. For instance, a configuration frequently observed on surfaces (e.g., in growth or ion 
bombardment) is one where the decaying island is on top of another larger island, forming 
a multilayer structure. Focusing on the decay of the top layer, the outer boundary of our 
treatment is the perimeter of the second layer. Because the energy barrier for interlayer 
transport of atoms is frequently high, the rate coefficient Kj, may be very small for this 
multilayer case. This situation, which we call the "outer-boundary limited" case, satisfies 
approximation (iii-a) (constant C). In this case the equilibrium density at the outer bound- 
ary is pr = as given by the Gibbs- Thomson relation. However, if the lower island is 
much larger than the upper one then the approximation (i)(-R ^ r) is also reasonable, so 
that pn ~ Poo and Eqs. (pTsD apply with C = KbR. Thus the outer-boundary limited case 
falls into the n = 2/3 mathematical limit. This configuration was considered by Villain,!^ 
who also found an exponent of n = 2/3. Note however that Villain's approximations lead to 
the diffusion limited case discussed above. Our outer-boundary hmited case gives the same 
exponent, but through a different mechanism. 




(15) 
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A final surface configuration of interest is tlie "island in a pit" case, where an island is 
within a monolayer pit. The outer boundary is well defined, similarly to the previous case. 
Here, however, the atomic processes at the outer boundary are the same as those taking place 
at the inner boundary, so that Kf, = Ki. There is no mathematical power law limit for this 
case, as neither approximation (iii-a) nor (iii-b) apply. We can, however, give a criterion as 
to whether the decay is dominated by the diffusion process or the interface transfer process. 
This criterion is similar to the cross over between the interface-transfer limited and diffusion 
limited cases above. This transition radius is defined by the condition that the second 
term of Eq. (|TTp is equal to the sum of the first and third, yielding the transcendental 
equation 

(16) 




K,R 



The left hand side of Eq. (|I6]) has a maximum value of approximately 0.278; if the right 
hand side is larger than this, then there is no transition radius and the entire decay is 
interface-transfer limited. Otherwise, the decay is interface-transfer limited for very large 
and very small islands, but diffusion limited between. The physical interpretation is that 
when the island is almost as large as the pit, diffusion of atoms across the narrow terrace 
between them occurs rapidly, leading to interface-transfer limited behavior for large islands. 
For very small islands, the short island perimeter leads to interface-transfer limited decay 
through lack of detachment sites. 

Finally, a comment on how experimentally controllable parameters should be expected 
to influence feature behavior. The overall rate of decay [see the prefactors of Eqs. ( p.3b|) 
and (|14b|) ] depends on many parameters, and therefore could vary widely from material to 



material. For any particular material, temperature will play a central role in controlling 
the decay rate, primarily through the temperature dependence of p^. This is because 
the equilibrium density results from adatoms being thermally activated to detach from step 
edges, and therefore will depend exponentially on inverse temperature. The effect adsorbates 
would have on a system is another particularly interesting influence on the decay rate.iiii 
Adsorbates lower the free energy of a surface. If they also lower the relative free energy of 
individual substrate adatoms, this could substantially increase p^, leading to an increase in 
the decay rate. Adsorbates would also be expected to lower the step free energy, improving 
the approximation (ii)(GT expansion) [see Eq. (0)]. However, if the adsorbates interact 
with the surface chemically through directional bonding, then such simple arguments are 
insufficient; other considerations, outside the scope of this paper, must be included. Finally, 
adsorbates could also affect rate constants by simply getting in the way of the substrate 
atoms, as mentioned above in connection with Eq. (p^Sj). 



IV. COMPUTER SIMULATIONS OF ISLAND DECAY 



The continuum model presented above is useful for identifying important factors con- 
trolling the stability and decay of surface features. However, it makes many simplifying 
assumptions, in particular about the proper macroscopic description of microscopic behav- 
ior. To approach the problem from another direction, we may examine island decay for 
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a generic system based on simple energetics, where different microscopic processes are ex- 
plicitly included. The results may then be compared to the thermodynamic model to see 
how well the latter captures the behavior of the system. To this end we have pursued com- 
puter simulations of island decay. One goal of these simulations is to verify that the above 
continuum treatment is applicable to nanoscale systems based on simple energetics, and in 
particular to test whether the limiting regimes of decay predicted by the thermodynamic 
model can be reproduced in simulations. We also hope to illuminate what microscopic 
mechanisms play the dominant role in the mass transfer process, and how they relate to the 
thermodynamic parameters. 

Our computer code simulates atoms moving on a square lattice using the solid-on-solid 
model, which does not allow overhangs. The system is rectangular with periodic boundary 
conditions. Atomic moves are classified according to the initial horizontal coordination of 
the moving atom (0 to 4 for the square lattice), the coordination the atom would have 
if the move were made, and whether the vertical motion is up, down, or absent. Taking 
into consideration that some moves are geometrically disallowed, this yields 48 types of 
moves. The rates for these move types are determined by an Arrhenius form, in which a 
single attempt frequency v and energy barriers for each type of move are parameters in the 
simulation. Note that the attempt frequency applies to each move separately, so that each 
atom attempts to move at the rate 4i/. In order to ensure that detailed balance is satisfied, 
these barriers are constrained such that the system obeys the bond counting Hamiltonian 

?i = JI^-B^ (coordination) , (17) 

sites 

where is the energy of a single bond. This constraint leads to 26 independent barriers, 
along with the bond energy and attempt frequency, as parameters defining the system. A 
dynamical Monte Carlo algorithmic!^ is used in which every computer iteration results in 
an atomic move, the simulated time being advanced probabilistically. This algorithm is very 
efficient relative to standard Monte Carlo methods, as no attempted moves are rejected. 

Although we are studying the generic behavior of surfaces, we wish to have barriers which 
are reasonable from a physical standpoint. Therefore, the energy barriers for in-plane moves 
used in the following simulations are based on barriers for the Cu(lOO) surface as calculated 
with the Finnis-Sinclair atom embedding model.il In Ref. ^ barriers were calculated for 
in-plane moves with all 2^*^ possible configurations of nearest and next-nearest neighbors 
to the initial and final sites of the moving atom. Exchange type moves were considered, 
but were found to be unimportant for this surface. For each combination of initial and 
final coordination we averaged an appropriate subset of these 2^'^ barriers to obtain the 
barriers for use in our simulation. This procedure clearly does not retain the details of the 
Cu(lOO) surface. However it does yield physically reasonable numbers, and in particular 
results in a bond energy consistent with the original barrier set. The energy barriers we use, 
corresponding to a bond energy of = 0.341 eV, are shown in Table |. 

The barriers for interlayer movement were set very high (^ 100 eV) in order to prevent 
the occurrence of such moves. In order for the decay to take place within a reasonable 
amount of computer time, the simulated temperature was a moderately high 1347 K, or 
60% of the critical temperature for this model. Note that the 2D critical temperature is a 



13 



TABLE I. Energy barriers for intra-layer atomic moves, in units of eV 



Coordination 




Final Coordination 




0-fold 


1-fold 


2-fold 


3-fold 


0-fold 


0.697 


0.479 


0.328 


0.166 


1-fold 


0.820 


0.624 


0.450 


0.275 


2-fold 


1.010 


0.791 


0.591 


0.377 


3-fold 


1.189 


0.957 


0.718 


0.462 



lower bound for the roughening temperature of the three-dimensional system.El The attempt 
frequency used was v = 10^^ s~^, which sets the overall time scale. 

Elsewhere we have shownlll that this model does exhibit a Gibbs-Thomson effect. That 
is, smaller islands are in equilibrium with higher adatom concentrations. There we show 
that the parameters of Eq. (|I]) can be analytically calculated0'@ for this simulation to be 
= 3.58 X 10~^ and 7 = 0.117eV/a. However, at the simulation temperature of 
1347 K, the adatom concentration of the vapor phase is high enough so that its equation 
of state deviates significantly from that of a noninteracting gas. Therefore, the proper 
Gibbs-Thomson relation is much more complex than Eq. (|l]). Nevertheless, we foundEll 
phenomenologically for this model at these temperatures that the GT vapor enhancement 
can be well described by including a correction factor g = 1.59 in the exponent of a modified 
Gibbs-Thomson equation 

p"=p-'''Kv§) ■ ^^^^ 

For our decay simulations, a 100 a x 100 a lattice was used. The initial configuration 
was a circular island, 15 a in radius, centered within a pit of radius 40 a [as in Fig. ||(a)]. 
As expected for a system exhibiting a Gibbs-Thomson effect, the island disappears over 
time as the adatoms move from the island to the pit edge. Figure H shows configurations 
typical of those observed during decay; the island remains relatively circular, and remains 
relatively well centered in the pit. In any individual simulation, the area of the island showed 
significant fluctuations during the decay (see Fig. ^ inset). To obtain an averaged behavior, 
ten simulations were run with identical initial conditions except for the seed for the random 
number generator. Each decay represents roughly 15 hours of CPU time on an IBM RS6000 
model 550 computer. The time origin of each run was offset so that the islands vanished at 
the same time, and the sizes of the islands were averaged at each sampled time. Very early 
times, during which the system was coming into steady state, were ignored. Figure ^ shows 
the averaged results. 

We wish to make quantitative comparisons between these results and the macroscopic 
theory of Sec. |I[ In particular, we are simulating the "island in a pit" case discussed in 
Sec. [illB| . As discussed there, the outer boundary rate coefficient is Kh = Ki. Because the 



geometry of the system is completely specified, it is unnecessary to make assumptions about 
screening lengths. The outer boundary position R is very close to constant, as the area of 
the island is small compared to the area of the pit. Finally, the equilibrium density at the 
outer edge is given by the GT relation as pr = p"'^. Note that the pit edge has a negative 
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radius of curvature, which changes the sign of the exponent in the GT relation. Thus the 
equation governing the decay of the island in this system is [see Eqs. (|T0|)] 



dr 
Itt 



Q / 1 



1 



1 1 

T ^ R 



-1 



91^ \ 
RkTj 



(19a) 



K 



D 



D 



ln{R/r) 



(19b) 



Here we have included the phenomenological correction factor g in the modified GT ex- 
pression. The definition of K^, is explicitly repeated to emphasize that it has significant 
dependence on r for this configuration. 

Note that neither the approximation (i)(-R ^ r) nor (ii)(GT expansion) applies to these 
simulations. Specifically, the island is too large relative to the pit for (i) to be satisfied, and 
the island is too small for (ii) to be satisfied. Therefore, the complete decay equation above 
must be used for a quantitative comparison of simulation and analytic theory. 



V. DETERMINATION OF MACROSCOPIC PARAMETERS FOR THE 

MICROSCOPIC SIMULATION 



In order to apply the macroscopic decay equation to the simulations, it is necessary to 
determine the correct values of the parameters D and Ki. These may be determined through 
proper measurements of the simulation behavior together with analytic thermodynamic 
calculations. Portions of the following treatment are achieved through the use of a mapping 
of the lattice gas to the Ising model, for which a high field (low concentration) expansion 
is used to obtain the equation of state. The interested reader can find a more complete 



description of this mapping and its use in Ref. 31. 




(b) t = 0.5 xlO-^ 



fc) t = 1.0 XlO"^! 



FIG. 2. Snapshots of a simulated island decay: the initial configuration (a) and later stages 
at (b) 45% and (c) 90% of the total decay time. White regions are one atom higher than black 
regions. 
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In the limit of vanishing concentration, the diffusion coefficient can be directly determined 
from the attempt frequency and the energy barrier for zero initial and final coordination: 



D{p^O) = aVe--^*/'^^ = 2.47 x 10^ a'^/s . (20) 

However, as noted above, the temperature of the simulation results in an adatom concen- 
tration which is high enough for the system to deviate significantly from ideal gas behavior. 
Therefore the diffusion coefficient was measured directly with "enforced gradient" simula- 
tions. Simulations were run in which a flat terrace populated by adatoms was divided into 
three sections. The densities in the left and right sections were held constant at values 
differing by 4 x 10~^a~^ (a small difference compared with = 3.58 x 10~^a~^). The 
resulting adatom current density was then measured and used to calculate a diffusion coeffi- 
cient through Pick's Law, J = —DVp where J is the adatom current density. Figure § shows 
the measured D versus the average of the densities in the left and right regions. As should be 
expected, an increasing density with an attractive interaction leads to a lower diffusivity.0 
Although the low densities make it difficult to obtain good statistics, the diffusivity is clearly 
reduced by approximately 25% from the dilute gas value. 

The second parameter which must be determined is Ki. As noted in Sec. 0, the as- 
sumption in the macroscopic theory that Ki is independent of r should be checked for any 



given system. Equation (||) shows that Ki ^ J'^. In Refs. ^ and ^ we show that J'^ may 
be measured from simulations through microscopic considerations, and that it changes very 
little with island size. That is taken as evidence that the assumption that Ki is independent 
of r is likely to be satisfied. However, our attempts at describing the simulated decay using 
Ki calculated from in this way have so far failed to agree quantitatively with the observed 
island decay, leading us to develop the methods described in the following paragraphs for 




t(10 \ec) 

FIG. 3. Simulated island decay using barriers based on Cu(OOl), showing the island area versus 
time. The figure shows the average area of ten simulations referenced to time of island disappear- 
ance. Error bars represent standard deviations of the mean. The inset shows two representative 
individual decays. 
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determining Ki from purely thermodynamic arguments. This procedure below was carried 
out for only one size island, under the assumption that is independent of island size. 

Our method of numerically extracting the interface transfer constant is closely related 
to the shape fluctuation analysis of Khare et a/.@ They study how the equilibrium fluctu- 
ations of the island shape depend on the interface transfer constant, the surface diffusion 
constant, and the edge diffusion constant; we will consider the equilibrium fluctuations of 
the island size. However, we stabilize the island against changes in size by placing it in a 
finite box with conserved order parameter, while Khare et al. stabilize the island size by ex- 
ternally adjusting the chemical potential under conditions of nonconserved order parameter. 
The equations in Ref. ^ therefore do not apply directly to our results. 

The parameter Ki may be determined through its relation to the mobility of an interface. 
Because a current of adatoms detaching from an interface causes the interface to move, Ki is 
closely related to the mobility for interface motion under the driving force of a vapor density 
imbalance. In principle, this mobility could be directly measured by putting an interface in 
interaction with a vapor of density lower than the equilibrium density. The interface mobility 
could then be determined from the rate at which the interface recedes. In practice for a 
computer simulation, however, it is very difficult to keep the vapor concentration constant 
near the interface. Detaching atoms continually increase the vapor density, and this source 
must be counterbalanced by an adatom sink. However, the required strength of this sink 
is not known a priori. Indeed, it is even difficult to distinguish whether a density increase 
near the interface is due to net detachment or is due to a fluctuation in the vapor density. 

The interface mobility is closely related to the fluctuations of the interface position in 
an equilibrium system. By considering an equilibrium system, the need for an adatom sink 
is eliminated. However, one must be careful because of the nature of the interface-vapor 
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• measured 
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FIG. 4. The dependence of the diffusion coefficient on adatom density in the simulation, ob- 
tained from direct measurements of current densities under enforced concentration gradients. Open 
circle shows the theoretical dilute limit value. The arrows indicate the concentration of vapor 
adatoms required for equilibrium with each of the pit (left arrow) and the island (right arrow) at 
the beginning of the decay simulation. 
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interaction. When an atom detaches from the interface, the local vapor density is enhanced 
momentarily. Because the vapor has a finite diffusivity, this enhancement will persist for 
some time, influencing the subsequent behavior of the interface. This memory effect means 
that the interface position does not perform a random walk, and the fluctuations of the 
interface position are thus slowed. The strength of the memory effect depends on the vapor 
diffusivity, and also on the system size because of the long range nature of the solution 
to the 2D diffusion equation. Note that these considerations do not influence the average 
vapor density near the interface. The memory effect is due to time correlations between 
fluctuations of the vapor density and fluctuations of the interface position. 

We have not obtained an analytic form for the interface fluctuations in the presence of 
such a vapor memory effect. However, it is clear that the memory can be removed by causing 
the adatom vapor to equilibrate quickly on the time scale of interface motion. In simulation, 
this can be achieved by randomly redistributing all monomer adatoms to other monomer 
sites between every atomic move. This explicitly does not affect the interface mobility, 
because all types of atomic moves which involve the interface are unaffected; the adatom 
redistribution is done in a way such that it will never cause the creation of new nearest 
neighbor bonds. By thus removing the vapor memory effect, the interface fluctuations may 
easily be analyzed to obtain the interface mobility, and hence the interface transfer constant 

In fact, the vapor memory is not completely removed by this artifice, because clusters in 
the vapor are stable. However, the density of clusters in the vapor is small in our simulations. 
Also, each cluster will have some mean time to dissociation, so that the memory introduced 
by clusters of a given size decays exponentially. For these reasons, the effect of clusters 
should be negligible. 

In detail then, we simulate an island in equilibrium with an adatom vapor in a closed, 
conserved mass system. Note that by simulating an island, rather than a straight interface, 
any angular dependencies are automatically averaged. Because mass is conserved, the island 
size is stabilized at some equilibrium size req. When the island fluctuates to be either larger 
or smaller than this, the resulting change in the surrounding adatom density tends to restore 
the island size. This can be understood in terms of a thermodynamic potential for the system 
which is a function of island radius, $(r). Reference ^ discusses this potential at length.0 
The equilibrium radius req is at the minimum of this potential. 

The fluctuations in the island size can be seen as a random walk of the island radius, 
characterized by some interface diffusivity Di. Note that the diffusivity of the interface Di 
is distinct from the adatom diffusivity D. In general the behavior of the interface would be 
affected by the adatom diffusivity D through the memory effect described above. For these 
simulations, however, monomer rearrangement is employed to achieve fast vapor equilibra- 
tion, so that the interface diffusivity obtained is a property of the interface alone. 

The random walk of the island radius is biased by the thermodynamic potential $(r), so 
that it tends to return towards rgq. The strength of this bias is given by the thermodynamic 
mobility of the interface -Bj, which is defined by the equation 

dr ^ d$ 
dt ' dr 
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-5, 



dp I dr 



(P - n + 



(21) 



eq 



In the second equality we have Taylor expanded d^/dr in p about the equilibrium density 
p^^, for reasons which will become clear below. Note that the interface mobility Bi and the 
interface diffusivity may be related through the generalized Einstein relation Di = kT Bi. 

To characterize the fluctuations of a random walk in a potential, we approximate the 
potential as quadratic in r: 



^(r) 



1 . 
-c r 

2 ^ 



eqy 



dr"^ 



(22) 



eq 



where c gives the curvature of the potential well. The value of c can be calculated from the 
known formic of $(r), but will not be required for the following analysis. Note that the 
Boltzman distribution in a quadratic potential is Gaussian, 



Probability(r) oc exp 



-$(r)' 
kT 



exp 



-c r — r, 



eq/ 



2kT 



(23) 



The Green's function for diffusive motion within a quadratic potential is therefore also 
Gaussian, with an exponentially relaxing mean and standard deviation: 



G{r,t]ro) 



1 



exp 



2(t2 



(24a) 



Tea + (ro 



-Bid 



(24b) 



El. 

B,c 



1 



-2BiCt 



(24c) 



where and a are the time dependent mean and standard deviation of the Gaussian form. 
Recall that the Green's function can be interpreted as the probability distribution which 
evolves from an initial delta function distribution at tq. It can therefore be used to find the 
time correlation function for the interface position r. 



(KO)r(t))-<=-^e 



D, 



-Bid 



(25) 



the prefactor of which may of course be simplified with the Einstein relation. Thus Bi is 



found by measuring the time correlation function of the island radius (using r = yArea/vr) 
and fitting it with Eq. as shown in Fig. ^. 

The interface transfer constant Ki is then obtained from Bi by comparing the ways they 
each describe interface motion. Ki describes how an interface moves due to a vapor density 
imbalance [see Eq. (H)], 



dr 2nr 

'dt~n' 



I, = 27irK,{p'^ - p{r)) 



(26) 
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In the Taylor expansion of Eq. (|2T|), the interface mobihty was also described in terms of 
densities, so that by comparing these equations we obtain 



^2$ 



--B. 

dr dp 



eq 




(27) 



cqy 



The derivative dr/dp is given by purely geometrical considerations. Note that it is not 
necessary to know c in order to obtain Ki with this method; the product BiC is obtained 
directly from the decay length of Eq. (^5]). As a consistency check, c was calculated from 
$ as known from the equation of state for the Ising modeLo It was then used to check the 
initial amplitude of the time correlation function. 



VI. QUANTITATIVE EVALUATION OF DECAY SIMULATIONS 

The above procedures yield values of = (2.47±0.19) x 10^'^ a/s and, at the equilibrium 
concentration for the initial island, D = (1.75 ±0.2) x lO^a^/s. Using Eq. (0), these values 
together with the average pit radius of R = 38.5 yield transition radii of r[ = R/4550 ~ 
8 X 10~^a and r[ = i?/1.0037 ~ 38 a. Thus the island, which is initially of size r = 15 a, 
remains well within the diffusion limited regime throughout the decay. Figure ^ shows the 
results of numerically integrating Eq. (|19a|) using these parameters. The agreement with 



the simulation results is very good. As expected based on the transition radii, the predicted 
decay is extremely insensitive to the value of the interface transfer constant. The deviation 
at very small islands may be due either to a breakdown in the phenomenological expression of 
the modified GT equation, or to a breakdown in the assumption that Ki is size independent. 
A breakdown of the quasi-steady-state assumption seems unlikely, since even at the highest 
recorded interface velocities each monomer makes over 1000 moves in the time that the 
interface recedes by one lattice spacing. 
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FIG. 5. The time correlation function for the interface position r, showing the function calcu- 
lated from one equilibrium simulation and a least squares fit of the form in Eq. ([25|). 
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In order to test the macroscopic theory in a regime where the interface transfer process 
is important, a second set of ten decay simulations were run in which the diffusion rate was 
greatly increased. This was achieved by reducing the barrier for moves from 0-fold to 0-fold 
coordination. That barrier was changed from 0.697 eV to 0.0 eV. This unphysical sounding 
"barrier" is possible because our simulation uses the barrier parameters only to determine 
rates. Note that at room temperature the same change in the rate of these moves could be 
achieved with a barrier reduction of only 0.14 eV. The large magnitude of the barrier change, 
required to obtain the following results, is necessary here because of the high temperature 
of the simulation. 

As determined by the methods described above, the diffusivity is even more concentration 
dependent than in the previous case. Figure ^ shows that at the vapor density which is in 
equilibrium with the initial island, the diffusivity D = 5.0 xlO^^ a^/s is one half of the dilute 
limit value of D = 1.0 xlO^^a^/s. As the island shrinks, the concentration increases and 
the diffusion coefficient decreases significantly. If we are to describe the entire decay with a 
single effective diffusion coefficient, it will be yet lower than 5.0 x 10^^ a^/s. 

Since non-diffusion barriers were unchanged, the same value of Ki as above applies here. 
Using these values in Eq. ([T6|), the left hand side is D/{KiR) = 0.317. As a result, Eq. (p!6D 
has no roots, and the interface transfer process is the most important process in limiting 
the decay rate during the entire island decay. Note however that this decay is not far into 
the interface-transfer limited regime, so that the diffusion coefficient is still expected to be 
of importance. 

Shown in Fig. ^ is the averaged decay from simulation, as well as analytic predictions. In 
this case the analytic results are sensitive to both D and Ki as expected. Note that the decay 
rate is two orders of magnitude larger than in the previous case. The relationship of interface 
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FIG. 6. The dependence of the diffusion coefficient on adatom density in the modified simu- 
lation which has a faster diffusion rate, obtained from direct measurements of current densities 
under enforced concentration gradients. Open circle shows the theoretical dilute limit value. The 
arrows indicate the concentration of vapor adatoms required for equilibrium with each of the pit 
(left arrow) and the island (right arrow) at the beginning of the decay simulation. 
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velocity and diffusivity is similar to the previous case, supporting the use of the quasi-steady- 
state assumption. The dashed curve, which uses the value D = 5.0 xlO^^ a^/s, qualitatively 
succeeds in showing this increase in decay rate, but deviates as the island becomes smaller. 
The entire decay can be matched quite well by using the value D = 3.0 xlO^^ a^/s, as shown 
by the solid curve. This can be partially justified by the fact that the smaller island results 
in higher adatom concentrations, and hence a lower diffusion coefficient. However, because 
of uncertainties in the calculations of D and Ki, it is difficult to compare simulation and 
analytic theory closely. 

In Sec. |T| we saw that in treatments of systems with no clear outer boundary, under 
certain approximations, the diffusion limited and interface-transfer limited regimes are char- 
acterized by power law decay of the island area, with respective exponents of 2/3 and 1. For 
these simulations, there is a definite outer boundary which is near the island and the various 
approximations are not valid, so that power law behavior cannot be achieved. However, 
at the end of the decay in Fig. |^ the decay rate increases markedly, whereas in Fig. |^ the 
increase in decay rate is less pronounced. In this sense, the diffusion limited decay is "less 
linear" than the decay dominated by the interface transfer process. This distinctive feature 
is correctly reflected by the analytic predictions. 

In principle, simulations might be run further into the interface-transfer limited regime 
by further increasing the difference between the energy barriers for interface transfer and 
for diffusion. However such an approach is impractical, as the program must spend most of 
its time processing diffusion events, and therefore decay occurs very slowly in real computer 
time. Other, more artiflcial means of increasing the effective diffusion rate, such as the 
monomer rearrangement scheme described in Sec. |V|, might circumvent this difficulty. 




FIG. 7. Simulated island decay with fast diffusion, showing the averaged island area versus 
time, referenced to time of disappearance. Error bars represent standard deviations of the mean. 
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VII. CONCLUSION 



Isolated features on a surface will tend to decay in size due to the Gibbs-Tliomson 
effect. This decay implies that there is mass flow in the adatom vapor surrounding the 
island, flowing from the island to the surrounding environment. There are three processes 
in this mass flow to be considered: interface transfer at the island edge, diffusion across 
the surface, and incorporation into the environment at some outer boundary. For various 
surface parameters and various surface conflgurations, each of these processes may be the 
rate limiting one, so that there are three general cases of the island decay. 

With certain approximations, the limiting regimes of these cases lead to a power law 
dependence of the island size on time. For the interface-transfer limited regime, this leads 
to the island area decaying as A oc (to — t). For the outer-boundary limited regime, this 
leads to the island area decaying as A (x {to — t)^/^. For the diffusion limited regime, no 
power law is strictly valid, but under a screening length hypothesis the A oc (to — t)'^^^ power 
law can be argued. 

The existence of these power laws may make it possible to experimentally distinguish 
between the regimes by measuring exponents for observed feature decay. Both the outer- 
boundary limited regime and the diffusion limited regime could probably be observed on 
physical surfaces. The interface transfer regime is unlikely to occur on clean metal surfaces, 
but may be important on other surfaces or on metal surfaces covered with adsorbates. It 
should be noted, however, that it can be quite difficult to distinguish between the two power 
laws based on a single decay. Data must either have very low uncertainty in measured island 
areas, or consist of many independent decays (as in Ref. in order to be a good test of 
the exponent. 

Computer simulations based on the solid-on-solid model and simple energetics exhibit 
the island decay behavior. Comparison of simulations and the analytic theory demonstrate 
that this theory correctly captures the essential elements of island decay. Although some 
limiting regimes of the theory cannot be fully achieved in simulation, the simulations match 
the theoretical predictions well for an important range of surface parameters. Note that this 
is the case in spite of the fact that the theory is thermodynamic and macroscopic in nature, 
while the simulation involves features of atomic scale. We conclude that even at atomic 
scales the analytic theory includes all the physics necessary to describe the behavior of a 
surface governed by simple energetics. 

The relationship between the parameters governing the microscopic dynamics of the 
system (e.g., the energy barriers for atomic moves) and the macroscopic dynamics parameters 
D and Ki is quite complex. We have shown methods of obtaining the macroscopic parameters 
from measurements of the simulation behavior. Further investigations of the relationship 
between microscopic and macroscopic parameters would be of interest, since in principle the 
macroscopic behavior is completely determined by the microscopic characteristics. 
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